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Solutions of a kinetic equation are investigated which describe, on a nonperturbative basis, the 
vacuum creation of quasiparticle electron-positron pairs due to a strong laser field. The dependence 
of the quasiparticle electron (positron) distribution function and the particle number density is 
explored in a wide range of the laser radiation parameters, i.e., the wavelength A and amplitude of 
electric field strength Eq. Three domains are found: the domain of vacuum polarization effects where 
the density of the e~e^ pairs is small (the "calm valley"), and two accumulation domains in which 
the production rate of the e~e^ pairs is strongly increased and the e~e'^ parr density can reach a 
significant value (the short wave length domain and the strong field one). In particular, the obtained 
results point to a complicated short-distance electromagnetic structure of the physical vacuum in the 
domain of short wavelengths A < Aacc = n/m. For moderately strong fields Eq < Ec ~ m^/e, the 
accumulation regime can be realized where a plasma with a high density of e'^ quasiparticles can 
be achieved. In this domain of the field strengths and in the whole investigated range of wavelengths, 
an observation of the dynamical Schwinger effect can be facilitated. 

PACS numbers: 42.55.Vc, 12.20.-m, 41.60.Cr, 42.55.-f 



I. INTRODUCTION 

In the present work we investigate the response of 
the physical vacuum (PV) of QED to the influence 
of a time dependent, strong periodic electric field 
("laser field") in a wide range of field parameters, 
i.e., the field strength amplitude and the wave- 
length A (wc use the natural units h = c = = ^, 
ks being the Boltzmann constant), 

0<Eo<Ec. = mVe, (1) 
ir/m = ttXc = -^acc < A < Ao, (2) 

where Ec is the so-called critical or Sauter-Schwinger 
field strength, and Ac = SttAc is the Compton 
wavelength of elementary leptons with mass m and 
electric charge e = |e| (we focus here on electrons 
and positrons); Aq ^ 1 /im ^. To characterize the 
PV response we use either the distribution func- 
tion or the number density of quasiparticle electron- 
positron pairs (EPPs) created from the PV. Below 
we select three characteristic domains in the plane 
spanned by Eq and A: the region of non-increasing 
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creation (dubbed " calm valley" ) and two boundary 
regions of accumulation, where the EPP density in- 
creases up to saturation (for sufficient duration of 
the field). The accumulation regions correspond to 
limiting high fields Eq ^ Ec (and A ;2> Ap) or lim- 
iting short wavelengths A '--^ Ac (and Eq < Ec). 
Of course, there is also a region where these ac- 
cumulation domains do overlap. Different mecha- 
nisms of vacuum creation act in these domains and 
provide the corresponding features of the PV re- 
sponse. The vacuum excitations in the calm valley 
can be described pcrturbatively while nonperturba- 
tive approaches arc necessary in both accumulation 
domains. In the following, wc consider the action of 
a periodical field. But the obtained results allow also 
to estimate the role of separate Fourier components 
of an external field with complicated spectral struc- 
ture (as, e.g., in the cases of vacuum EPP creation 
by scattering of charged particles at high energies [l| 
and the dynamical Casimir effect Q). 

The methodical basis of our analysis is a system 
of kinetic equations (KEs) which is a nonperturba- 
tive consequence of QED in the presence of a spa- 
tially homogeneous, time dependent electric field [3- 
Q . These KEs arc intended for the description of the 
dynamical Schwinger effect of vacuum particle cre- 
ation in terms of quasiparticle vacuum excitations, 
having the quantum numbers of the electron and 
positron in the presence of an external field (mass, 
charge, quasi- momentum and quasi-energy), which 
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are called below as quasiparticles. After the ceas- 
ing of the laser pulse some of these quasiparticles 
become real electrons and positrons with their mo- 
mentum and energy lying on the mass shell. 

Thus, in the present work, vacuum polarization 
effects and vacuum particle creation in the time- 
dependent electric field of a standing wave formed 
in the focal spot of counter-propagating strong laser 
fields are analyzed in the framework of the quasi- 
particle representation. Such a quasiparticle pic- 
ture corresponds to the level of description with real 
EPPs in the out-channel. On the other hand, the in- 
formation about quasiparticle distributions is useful 
for different estimates of the secondary observable 
effects such as the generation of pair annihilation 
photons the birefringence effect (e.g., [l2|) 

and so on. The additional investigation of the resid- 
ual EPP surviving after the ceasing of the laser pulse 
will be given elsewhere. 

The present work is organized as follows. In 
Sect. In] the basic KE and the set of limitations for 
their applicability are presented. The results of nu- 
merical investigations of the vacuum EPP creation 
kinetics for the case of a linearly polarized laser field 
are discussed in Sect. Illll for different domains of the 
external field parameters in the ranges ([T|) and 
In the calm valley of the laser radiation parameters, 
the same results can be reproduced analytically in 
the low density approximation by using the peculiar 
perturbation theory (see Sect. Iivp . The summary 
of the results is given and discussed in Sect. |Vl Here 
the actuality of development of a nonperturbative ki- 
netic theory of vacuum EPP creation for the case of 
non-quasiclassical external electromagnetic fields is 
emphasized for cases where the quantum field fluctu- 
ations are essential. It is also shown that both accu- 
mulation domains are promising for an experimental 
observation of the dynamical Schwinger effect. 



II. BASIC KINETIC EQUATIONS 

It is well known that vacuum particle creation 
in strong electromagnetic fields is possible if one 
or both of the field invariants Ii — — and 
I2 ~ EH are no n- vanishing [l3-[l7|. Such condi- 
tions can be realized, e.g., in the focal spot of two 



or more counter-propagating laser beams, where the 
electric field is spatially homogeneous over distances 
~ A and time dependent. For very strong fields 
Eq Ec one can then ex pec t the creation of a real 
electron-positron plasma [17h2]| out of the clcctro- 
magnetically polarized PV. In the subcritical field 
regime Eq <^ Ec the question is mainly about the 
creation of a short-lived quasiparticle EPP that ex- 
ists during the course of the field action (the dy- 
namical Schwinger effect) [l,!!^. The possibilities 
for observing this kind of PV response with the help 
of secondary effects such as the radiation of anni- 
hilation photons are investigated intensively at the 



present time. 

In this work, we study the PV response to a 
periodic model laser pulse with linear polarization 
A'"(t) = (0,0,0, A(i)) (Hamiltonian gauge), where 

A{t) = {Eo/iy)cos{iyt), E{t) = Easm{iyt) . (3) 

We are going to investigate the PV response as a 
function of the angular frequency ly ~ 2tt/X (or the 
wavelength A) and the amplitude Eq. For this aim 
we use the exact nonperturbative KE of the non- 
Markovian type for the one-body EPP phase-space 
distribution function f{p,t) obtained in Q, 

/(p,t)-A(p,t) J dt'A{p,t')[^-f{p,t')]cos29{t,t') 

(4) 



where 



A{p,t) = eE{t)e^/uj^{p,t) 



(5) 



is the amplitude of the EPP excitations with the 



quasi-energy w(p,<) = ^ e\{p) -f {p\\ — eA{t)Y and 

£_L = (w^ + p1)^^^ as the transverse energy; the 
quantity 



e{t,t') 



dr w(p,t) 



(6) 



is the high frequency phase. The distribution func- 
tion in the quasiparticle representation is defined 
relative to the in- vacuum state, f{p,t) = {in \ 
a^{p,t)a{p,t) I in), where a(p,<) and a~^{p,t) are 
the annihilation and creation operators in the quasi- 
particle representation. The existence of the quasi- 
cncrgy oj{p,t) suggests that the quasiparticle exci- 
tations arc not on the mass shell. For the general- 
ization of the KE ([4|_to arbitrary electric field po- 
larization, see Refsni,!!!,!!!. 

The KE (HI is equivalent to the system of ordinary 
differential equations 



f = -Aw, 

u = A(l - 2/) - 2ujv, 
V — 2u)u, 



(7) 



which is convenient for the subsequent numerical in- 
vestigations. The functions u{p,t), v{p,t) describe 
the vacuum polarization effects. Initial conditions 
at t = are f = u = v = 0. 

It is assumed that the laser electric field E{t) ([3]) is 
quasiclassical. This means that the photon number 
with the frequency v = 2Tr/X must be rather large in 
a volume of the order A"^ . This condition is fulfilled 
m,!!! for 



Eo > (i^/27r)2 



(8) 



i.e., in the quasiclassical (QC) domain. In the quan- 
tum (Q) domain for Eq < (z^/27r)^ it is necessary to 



take into account the quantum fluctuations of the 
external electromagnetic field and a corresponding 
generalization of the KE Q is required. 

Below it will be shown that the features in the 
behavior of the PV response for the weak field case 
Eq <^ Ec become apparent just in the Q domain, 
where the applicability of KE (|4|) breaks down. An 
exception is the strong field domain Eq ^ Ec, where 
the inequality ^ is valid. However, in spite of this, 
we will investigate its solutions here assuming that 
the external field can be treated as some quasiclassi- 
cal background. Thus, these extrapolated solutions 
have the character of a preliminary forecast. 

The adiabaticity parameter |27l - [29| 

is introduced for separating the domains of influ- 
ence of two mechanisms of vacuum particle creation: 
tunneling for 7^1 and multiphoton for 7 3> 1 pro- 
cesses. The point 7 = 1 separates these domains; it 
corresponds to the boundary curve EcXc ~ EqX- 

It is assumed that the field ([S]) is switched on at 
to = (we ignore here the error brought in by such 
an instantaneous switching-on [TtI . |30|V This way 
of switching on the external periodical field is rather 
standard in the framework of the discussed problem 
pTl ITsl . Below we will consider also the method 
of adiabatic field switching on. Generally speaking, 
the obtained results depend on the way of the field 
switching on. This is due to the occurrence of high 
harmonics in the solution spectrum of the KE 
As a rule, the high harmonics, which accompany 
the fast (non-adiabatic) switching on of a laser field, 
have no essential influence on the solution spectrum. 

For a comparative investigation of EPP produc- 
tion as a function of the flcld parameters A and 
Eq it is convenient to use some indicator charac- 
teristics. For this purpose we introduce here the 
maximal EPP number density rimax — nit — T/4) 
which is related to the maximal distribution func- 
tion /max(p) = f{p,t = T/4) obtained from solving 
the KE (HI) by integration over the momentum space 

n{t)^gl ^fip,t) , (10) 

where g = Ais the degeneracy factor, taking into ac- 
count the spin and charge degrees of freedom. This 
definition is used below in the case of absence of 
accumulation effects for which the number density 
EPP grows rapidly. In Fig. [1] the above discussion is 
summarized in the "landscape" of the laser param- 
eters Eq and A. 

III. SHORT-DISTANCE 
ELECTROMAGNETIC STRUCTURE OF 
THE PV 

A numerical investigation of the KE (|3|) demon- 
strates the complicated behavior of the distribution 
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FIG. 1: (Color online) Electron-positron pair production 
in the "landscape" of field strength Eo and wavelength 
A of laser collider experiments, with a few examples 
for comparison (XFEL [11], Astra HI], ELI [H). The 
hatched areas show the regions of accumulation due to 
the strong field (F) or resonance (R) mechanisms, resp. 
The black lines (dotted, dash-dotted, dash-double dot- 
ted) show three examples for peak pair densities rimax. 
The violet dashed line depicts the limit between the 
quantum (Q) and the quasiclassical (QC) domain which 
itself is limited by the red solid line (adiabaticity index 
7 = 1). 



function which shows a series of qualitative modifi- 
cations where the laser radiation parameters Eq and 
A are varied in the boundaries ^ and It is 

important to note that the character of the distri- 
bution function evolution depends strongly on the 
selection of the momentum representation: /(p, t) 
accomplishes some oscillations swinging along the py 
axis ^ (the direction of the field ([3])) as a whole with 
the amplitude I/7 simultaneously with the "breath- 
ing" mode, when its amplitude and form are altered. 
The transition to the kinematic momentum elimi- 
nates these oscillations and keeps the "breathing" 
mode only. 

The characteristic domains of dynamical behavior 
of the distribution function f{p,t) and the number 
density of EPPs (fTO|l in the full range of the laser 
radiation parameters ([1]) and ([2]) are shown in Fig. [1] 
to be discussed step by step in what follows. The re- 
gion outside the hatched accumulation domains due 
to strong field (F) and resonance (R) mechanisms of 
pair production is the calm valley. The boundary 
between quasiclassical (QC) and quantum (Q) do- 
mains of the electric field Q is given by Eq = 1/A^ 
and the line 7 = 1 separates the multiphoton and 
tunneling domains. Surprisingly, rimax is A indepen- 
dent practically in the whole calm valley, i.e. outside 
the accumulation domains R and F. Its dependence 



^ This follows from the definition of the kinematic momen- 
tum P = p„ — eA{t) and the construction ((Sjl of the laser 
field. 
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on the strength of electric field Eq is shown in Fig. [21 
The symbols on Fig. [1] depict also the basic param- 
eters in the focal spot of the existing [sil, [32| and 
planned jssl - lssj laser systems as, e.g., Astra, XFEL 
and ELI. 
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FIG. 2: The dependence of rimax on the field strength 
Eo in the calm valley. 

In the following subsections we discuss the three 
characteristic domains of Fig. [1] separately. 

A. The calm valley domain 

This domain corresponds to the most simple be- 
haviour of the distribution function characterized by 
the absence of an appreciable production rate for 
EPP when averaged over the period of field oscilla- 
tion. 

It is interesting to consider the initial behaviour of 
the distribution function /(p) immediately after the 
field is switched on. In Fig.[3]wc present the typical 
form of /(p) at the time t = T/A, where T ~ 2tt/i/ 
is the period of the laser field oscillation. 

This distribution is anisotropic and can be approx- 
imated by a two-temperature Boltzmann distribu- 
tion 

/eq(pi,P||) ~exp{-ei/Ti -£||/r||} , (11) 

where T±_ is the transversal and Ty the longitudinal 
temperature while e||(p) = y^m^ + p| is the longi- 
tudinal energy. Fig. [4] demonstrates the similarity of 
the distributions /(p) and feq{p) (fTT|) and the rela- 
tion for the temperatures ^ = T\\/T^. However, a 
more satisfactory approximation for the pair distri- 
bution function is 

/in(p)^ei(p^)/c.o'(p) . (12) 

The difference between the distributions (fTT|) and 
([T^ appears in the high-energy tails. The estimate 
((T2|) can be obtained on the basis of the KE ([4]) in 
the low density approximation in the multiphoton 
domain 7 3> 1 (see also Sect. lIVp . In the framework 




FIG. 3: The EPP distribution function /(px.py) for the 
laser radiation with A = 1 nm and E = 10'* E^ at the 
moment t = r/4. 



of the approximation ()lip the temperature param- 
eters Tl and Ty are universal and do not depend 
on the laser field parameters A and Eq. In these 
assumptions the distributions pT|) and p2|) are in- 
dependent of the wave length A and proportional to 
i?g. Moreover, in the initial stage these estimates re- 
main valid regardless of the form of the electric field 
pulse. This holds in particular for the Sauter poten- 
tial, where the solution is well known [itI l40ll4l|. 

The excitation mechanism of the anisotropic spec- 
trum of EPPs in the early stage of the evolution is 
fully determined by the dynamics of vacuum cre- 
ation and, specifically, by the structure of the am- 
plitude ([S]). There the anisotropy effects are repre- 
sented by both, the transverse energy £_l(pj_) and 
the quasienergy aj(p, t). The same anisotropy factor 
^ = 2/3 can be obtained on the basis of the distri- 
bution ((T^ . when the relation of derivatives of the 
projections fin{p\\,P± = 0) and /j„(p|[ = 0,p_l) with 
respect to ey and ej_ is considered. 

Thus, in the case of a short pulse, the anisotropic 
distributions (|lip and will lead to an elliptic 
flow of EPPs, which is compressed in the direction 
of the electric field. 

The estimates (fTT|) and p2|) are modified when 
other mechanisms of vacuum creation of fcrmions are 
realized (see, e.g., [1^). This can lead to a change 
of the anisotropy factor ^. Considering, instead of 
fcrmions, the case of bosons, the situation changes 
radically: the equilibrium-like distributions of the 
type (fTTj) and (fT2l) a re replaced by the ones for strong 
nonequilibrium 137|, see Fig. [5l 

Thus, one can say that the initial distribution of 
the particles created from vacuum is well defined in 
the general case. In particular, the conclusion about 
the existence of equilibrium-like universal distribu- 
tions can become important for the theory of quark- 
gluon plasma generation in heavy-ion collisions. 
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FIG. 4: Fitting the KE solution /(p, t = T/A) of Fig. H 
(dash-dotted line) with an anisotropic equilibrium Boltz- 
mann distribution with T± = 0.3 m (upper panel) 
and = 0.19 m (lower panel) as well as with a power- 
lilce distribution (|12|) . 




FIG. 5: Momentum space distribution of bosons for the 
same conditions as in Fig. [S] 



A few remarks are in order. Apparently, the tem- 
perature parameters Tl and T|| can not be inter- 
preted as the Unruh temperature Tu = a/27r [isj . 
since the acceleration a ~ —{e/m)A{t) ^ Eq would 
imply a field dependence of T± and Ty which is not 
observed in our numerical investigations. Finally, 
some anisotropic initial distributions of particles cre- 
ated from vacuum under the action of an electric 
field have been introduced long ago (see, e.g., j36j). 
but they have artificial character. 

Such a simplifying picture of the phenomena in 
the calm valley is valid approximately in the ini- 
tial stage of the process for t < T. The high fre- 
quency harmonics become essential at t ^ T and 
the structure of the distribution function gets very 
complex. However, the amplitude estimate remains 
valid and the EPF production rate averaged over 
a period (i.e., the pair creation per unit time) is 
absent, < / >t~ 0. The dependence nmaxC-Eo) as 
depicted in Fig. [2] reveals that in the whole calm val- 
ley domain holds n.max(£'o) ~ Eq up to the strong 
field accumulation region (the shaded area "F" in 
that figure) where the numerical calculation is com- 
plicated. The distribution function f{p,t) shows a 
breathing oscillation with the amplitude /max and 
a frequency being twice that of the laser field (see 
Fig. [6l dashed line) . 



n(t) [fm-^] 
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t[T] 

FIG. 6: Time dependence of density n{t) in the weak 
field regime evidencing the transition from the oscillating 
mode for A = 0.01 nm to the linearly growing (in the 
mean) mode at A = 0.001 nm. The latter behaviour is 
termed here as "accumulation". 

These features of the distribution function be- 
havior are reproduced very well analytically on the 
quasiparticle level in Sect. IIVI The residual number 
density of EPF is very small in comparison with the 
breathing mode. 

The calm valley is bounded from above and to 
the left by the two accumulation domains F and R 
(see Fig.[T]), where the amplitude of the distribution 
function increases with the lapse of time and the av- 
eraged EPF production rate < / >t becomes appro- 
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ciable. Apparently, a glimpse of accumulation effect 
is presented in the calm valley too. However, this 
effect is negligibly small for i/ ^ m and Eo <^ Ec- 
Thus, the accumulation effect becomes dominant in 
the F and R domains only. 



B. The resonance accumulation domain 

This domain starts in the short wavelength do- 
main from the point Aacc ~ ttAc — 0.0012 nm, cor- 
responding to the energy v — 2m, which peaks out 
from the monochromatic external field for creation of 
an EPP. Wc will denote this process as one-photon 
pair creation. Below we will use this term in the 
analysis of solutions of the KE (|4|) , for which the spa- 
tial inhomogeneity effects are negligible. Here, in the 
point Aacc the accumulation mechanism is switched 
on sharply. As an example, the initial growth of the 
density in the course of time is depicted in Fig. [6] 
for A = 0.001 nm and Eq = lO"'* E^. At the later 
stages the degeneration effect developes. Here the 
distribution function reaches its maximal value and 
after that it performs oscillations which are damped 
asymptotically (see Fig. [7]). In addition, the 2i/- 
dependence of the breathing mode is conserved in 
the R domain. 
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FIG. 7: The limitation of accumulation due to the sat- 
uration of /(p) at large times for A = 0.001 nm. The 
saturation is achieved faster for stronger fields. 

The shape of the distribution function is changes 
rapidly. Later on, for A < Aacc, it leads to a collapse 
of the distribution function in a thin spherical layer 
(EPP bubble, see Figs. [8] and [9l). The thickness 
Ap and the size of this layer depend on the field pa- 
rameters A and Eq- For example, for A = 0.001 nm 
and Eq = 0.24 Ec we have Ap ^ 0.1 m. With grow- 
ing field strength at a fixed frequency the maximum 
value /max(p) increases in the degeneration condi- 
tion, /max(p) — 1- The occupation of the EPP bub- 
ble increases also (Fig. ^ under this condition (ap- 
pearance of some thin wave structure arises within 
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FIG. 8: The distribution function f{p±,P\\) in the R do- 
main at time t^T. Most of pairs are seen within a thin 
spherical layer centered at the origin of the coordinates. 




FIG. 9; The /(0,P||) dependence in the R domain at 
large times after saturation for different field strengths. 



the bubble). On the other hand, the decrease of the 
wavelength in the domain A < Aacc does not change 
the distribution picture essentially (see Fig. ITU)) . The 
radius of the EPP bubble in momentum space is 
equal to 



Ph 



ub = V {t^I^Y - ™^ for A < Aa 



(13) 



For V — > 2m, this radius goes to zero, Pbub -> (it 
corresponds to the condition v = 2m or the point 
A = Aacc)- Decreasing the wavelength is accom- 
plished by the growth of the EPP bubble size (see 
Fig. nop . The maximal value of the number density 
of the vacuum bubble in the saturation state (see 
Fig. [TT|) and the time period necessary for achieving 
the saturation (see Fig. [T2]) show that the first pas- 
sage time of the saturation grows indefinitely with 
decreasing A < Aacc- From Fig. [T] it follows that the 
point Aacc lies on the left of the critical line 7 = 1 
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FIG. 10: The /(0,P||) dependence in the R domain at 
large times after saturation for different wavelengths. 
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FIG. 11: The dependence of Umax on wave length A in 
the R domain for Eq = 0.1 Ec. 
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FIG. 12: The dependence of the saturation time on 
the wavelength A in the R domain for Eq = 0.1 Ec 



in the domain 7^1, where the multiphoton mech- 



anism of EPP excitation is operative. Apparently, 
one can expect that the multi-photon process can 
be changed to a few-photon one in the lower left 
part of this figure. As the result, the few-photon 
domain is limited by the two- and one-photon pro- 
cesses. For weak fields E <C E,., the one-photon 
process of the EPP creation dominates (this con- 
clusion is confirmed by the analytical calculations in 
Scct. lIVp . The two-photon pair creation process (the 
inverse of the Breit- Wheeler process [Isl, HB, [11] ) is 
switched on for rather strong field Eq ~ Ec begin- 
ning with X2-y = Ac = 2Aacc, see Fig. [T31 This 
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FIG. 13: Transient region from calm valley to F accu- 
mulation: the characteristic for F domain with linear 
growth of mean density is here caused by the subcritical 
field Eo = 0.5Ec, at the wavelength for the opening of 
the two- photon production channel A = Ac = 2Aacc. 

domain is characterized by the fast growth of the 
EPP density and it is contained in the F accumu- 
lation domain. Surprisingly, the KE ([4]) is sensitive 
to this effect which is valid in the case of a quan- 
tized electromagnetic field. Let us remark that the 
one-photon mechanism of EPP excitation prolongs 
to act in all domains of short wavelengths A < Aacc- 
After having discussed implications of the purely 
periodic field ([3]) we turn now to two other exam- 
ples of the time dependence. The initial condition is 
fit ^ -oo) = 0. 



1. Gaussian envelope 

Let us investigate first the PV response to the 
action of a pulsed field with Gaussian temporal en- 
velope 



E{t) = Eq sin {i^t) exp {-t^/r'^) 



(14) 



The comparison of the accumulation effects for the 
periodical field ([3]) with A = 0.001 nm (dashed line) 
and the corresponding field pulse (|14p with tq = 10 
T is shown in Fig. [M) It allows to draw the conclu- 
sion that the form of the field pulse is rather essen- 
tial but it does not change qualitatively the picture 
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of the effect. The switching-off process of the pulse 
(jl4p is accompanied by a stepwise reduction of the 
EPP density down to some residual level, see Fig.[T4l 
which should be compared with Fig. [51 Let us re- 
mark that in the regime of saturation the residual 
EPP density can surpass considerably the amplitude 
of the breathing oscillations while in the calm valley 
the situation is opposite. 




FIG. 14: Transition from the oscillating mode (for A = 5 
pm) to the accumulation mode (for A = 1 pm) for the 
pulse ([Till with TO = 10. 



2. Sauter pulse 
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4x10 



2x10 




FIG. 15: The residual pair density function of 

Ts for the pulsed field (|16p . 
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Now we consider action of a smooth pulse field 
with the Sauter potential 

A{t) = -EoTs tanh {t/rs) , (15) 
E{t) = Eo/ cosh^it /ts) . (16) 

Similar to the field (fT4|) one can clearly define here 
the residual density riout = "^1(1 — oo). The typical 
resonance picture is revealed here with the maxi- 
mum at Ts ~ 1/2™, see Figs. [TSl and [TBI For 
Ts ^ l/2m a sharp decrease of the PV response 
takes place and the vacuum EPP creation practically 
ceases. This result corresponds to the conclusion of 
Ref . . The characteristic value = 1/ 2m allows 
to speak here about the infiuence of one-photon pair 
production process on the level of the kinetic de- 
scription. This feature in the solution of the KE (|4]) 
with the potential ([T5|) can be confirmed also by an 
analysis of the exact solution of the problem [4l| . 

The accumulation effect in the R domain can be 
explained if one takes into account the presence in 
the solution of the KE of a parametric resonance 
in the neighborhood ~ 2m and of all combina- 
tional frequencies ni> and 2ml, where n and I are 
the integers and the factor 2 is stipulated by the 
structure of the high frequency phase ([S]). This is 
compatible with a parity odd distribution function 



FIG. 16: The time dependence of pair number density 
for the pulsed field (|16p for three values of the pulse 
length parameter rg. 



f{p,t) under time reversal. It leads to the appear- 
ance of contributions of the type 

— sm(v — 2m)t ■ sin (z^ -|- 2m)< . (17) 

— 2m 

For v 2m a secular term appears here that results 
in a linear growth of the distribution function. The 
subsequent evolution is accompanied by a growth of 
the EPP number density and saturation as a result 
of the action of the statistical factor in the KE 
Just this picture is observed for the numerical solu- 
tion in the R-accumulation domain. The given inter- 
pretation was generated by the perturbation theory 
in the low density approximation, see Sect. IIVI 

Let us remark also that the motion along the A- 
axis on the side of short waves is accompanied by a 
sequential replacement of the vacuum creation mech- 
anisms: the tunneling mechanism acts in the limit 
of stationarity (A oo) and slowly alternating field 
and requires an infinite photon number from the ex- 
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tcrnal field reservoir. This mechanism is then re- 
placed by the multiphoton (7^1) and few-photon 
ones and finally turns into the one- and two-photon 
pair production. It is necessary to underline that 
the conception of the "photon" is considered here 
in the framework of the accepted model of the spa- 
tially homogeneous electric field acting in the focal 
spot of counter propagating laser fields with the lin- 
ear size ^ A that allows to use in the description of 
the absorption processes the energetic condition only 
(e.g., V = 2m for the one-photon e~e+ production 
process). 




0.61 



^^^nr^o -30 -40 -bO -et 

P„ [m] 



C. The strong field accumulation domain 



After the interlude on finite pulses let us return 
to the discussion of further effects related to the 
periodic field ([3]). For the first time, the accumu- 
lation effect in a strong laser field (domain F in 
Fig. [T|) was discovered in Ref. [4^ on the basis of 
a numerical solution of the KE (|4]). The accumu- 
lation domain begins here with a rather high field 
strength Eq w 0.2Ec. The smooth breathing mode 
with a smooth phase space distribution (Fig. [3]) is 
replaced here by a strongly fragmented one with 
pronounced structures (Fig. [T7|). Apparently, in the 
long wavelength limit the usual accumulation mech- 
anism acts here. For the transition to the observ- 
able domain, the electron and positron must gather 
within a spatial extent of one Compton wavelength, 
an energy which is comparable to the energy gap, i.e. 
eEcXc = 2m. In the multi- and few-photon domains 
the other mechanism is switched on where the EPF 
creation is a result of the simultaneous absorption of 
n identical photons from the "photon reservoir" of 
the laser field. The probability for such kind of pro- 
cess, in principle, can be taken into account on the 
basis of an analogy with the theory of multiphoton 
ionization of atoms given, e.g., in Ref. [29| . 



The intersection region of the two accumulation 
domains R and F (the left top corner on Fig. [1]) is 
especially interesting. Here the nonlinear superpo- 
sition of two accumulation effects can lead to some 
new features in the vacuum EPF production. It is 
important that the applicability of the KE ([4]) is not 
violated in this domain. 

Thus, the accumulation processes in the R and F 
domains are the processes of the EPF excitation up 
to the saturation state and next increasing of EPF 
density in this state. Let us note that the vacuum re- 
sponse reveals the symmetry of the states in a weak 
{Eq -C Ec) and a strong {Eq ^ Ec) field, respec- 
tively: under the substitution / = 1 — / in the KE 
dH) it does not change its form and the case / <C 1 
is similar to the case / ^ 1. 



f(p,T) 
0.008 



0.006 



0.004 



0.002 




-29-5 .30.0- 
P,, [m] 



30.5 



31.0 



FIG. 17: The distribution function in the F accumula- 
tion domain at t = T (Eq = 0.5 Ec, A — 0.15 nm). 
Upper panel: general view; lower panel: detailed view 
highlighting the concentric dense layers. 



IV. LOW DENSITY APPROXIMATION 

The numerical results discussed in Sect.|lll]can be 
confirmed analytically in the multiphoton domain 
(7 ^ 1) of the calm valley. From the KE ([4| it 
follows in the low densit y ap proximation / <C 1 cor- 
responding to Eq ^ Ec [201 that 



./(P,t) = 



dt'X{p,t')e 



(18) 



This formula solves correctly the initial value prob- 
lem in the KE ([4]) providing a vanishing distribution 
function both at < — > to and £'0 ^ 0. In the leading 
approximation uj{p,t) — > aJo(p) = ^/rrl?+p^ , cor- 
responding to 7 3> 1, after the time integration, one 
arrives at 



/(p,0 = /(0)(p)+2/(2)(P)[l 



COS 2iyt\ 



(19) 
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where (a = e^/47r) 

/(0)(P) = T^aEl 
/(2)(p) = ^aEl 



V 



2^2 



2 



£ I 



(20) 
(21) 



In these calculations, the assumption of adiabati- 
cally switching-on the electric fields (1) at <o ^ ^oo 
was made. 

The leading role of the second harmonics (fT9|) was 
detected in the numerical solutions of the KE ([4]) 
(see Sect. imp and used in 0,[22| with the amplitudes 
/(o) and /(2) as the results of the numerical analysis. 

The Eq dependence of /max and nmax (see Fig. [2|) 
follows immediately from Eqs. ([^0]) and (PT]) . The 
time dependent component (|2ip does not depend on 
1/ for <C 771 (see Fig. [Tl calm valley); the con- 
stant component (j20p brings some small dependence 
in this domain. The calculations close to the poles 
in Eqs. ([20)1 and (|2ip require a more careful analy- 
sis. The Eqs. (PU)) and (PT|) indicate the one-photon 
mechanism of the EPF bubble generation in the 
resonance domain v ~ 2a;o(p) (sec Sect. IIIII Sub- 
sect. IIII Bp . In the considered approximation this 
bubble has a spherical form. 

The estimates of the density ([T0|l on the basis of 
Eqs. (201) and (EH) lead to the result 

n{t) = 77(0) + ^77(2) [1 - cos(2i/t)] . (22) 

In the leading approximation with respect to small 
v/uiQ it follows from Eqs. and ([21]) that 



1(0) 



77(2) 



I- i - 

A m \ lb J \m 

2 



3 a f Eq 



2 777 V 4 



(23) 
(24) 



The ratio of these expressions is equal to 77(o)/77.(2) = 
(5/72)(i//777)^ , which is close to the estimate given 
in 0. The frequency independence and Eq depen- 
dence of the amplitude n(2) of the number density 
oscillations agrees with the numerical calculations 
exhibited in Fig. [TJ The constant (residual) com- 
ponent density ([23]) was not observed previously in 
numerical calculations. The numerical integration 
shows that the estimates (|23p and ((24|) remain valid 
up to the quasiclassical boundary ly < 2777. The 
asymptotics at i/ — t- in Eqs. ([20]) and ([23]) leads to 
a vanishingly small residual EPF density. 

Thus, in the multiphoton domain of the calm val- 
ley the averaged production rate < / >t of EFF 
is negligibly small. Instead, there are vacuum oscil- 
lations induced by the external field and the small 
constant (residual) component. In principle, these 
vacuum effects can lead to some observable effects. 
For instance, they influence on the propagation of an 
electromagnetic probe signal. This can result in an 



anomalous absorption of this signal at the frequency 

27/. 

The perturbation theory used here in the quasi- 
particle representation is valid in the multiphoton 
domain 7 3> 1, when the swing amplitude (Sect. UlI)) 
is small. However, the numerical solutions of the KE 
([4]) show that these results are valid also in some 
parts of the tunneling mechanism region. At first 
glance, this fact is unexpected. Indeed, in this re- 
gion for subcritical fields Eq ^ Ec it is very well 
known that the EFF production rate is exponen- 
tially small [13, [3 ■ But this effect is not accessible 
to perturbation theory. On the other hand, the do- 
main of the long wavelength limit A > Aq is difficult 
for the numerical analysis and it was not studied in 
the present work. It is possible, that some special so- 
lution of KE ([4]) is present here which is compatible 
with the well known estimates fl3, UM ■ 



V. SUMMARY 

In the present work we have described the FV 
response to a monochromatic laser radiation in the 
landscape of the laser parameters ([T]) and ([2]). Three 
specific domains were observed here 

- the domain of the vacuum oscillations induced 
by the external field (calm valley), where the 
EPF production rate is very small, < / >t~ 
0, and the breathing modes are acting only, 

- the domain of the accumulation effect due to 
strong field (F), and 

- the resonant (R) accumulation domain, where 
the EPF production rate can reach significant 
values. 

Thus, it was shown that the behavior of the FV 
in the calm valley up to a certain critical point 
Aacc = ttAc is stable and does not depend on the 
wavelength A of the external electromagnetic field. 
This domain accommodates different mechanisms of 
vacuum decay into quasiparticle EPF excitations 

- tunneling (7 <^ 1, when the photon number of 
the external field is infinity or very large), 

- multiphoton processes (7 > 1) and 

- few-photon processes. 

The accumulation effect depicted in Fig. [l] ap- 
pears just cither in the multiphoton domain as a 
parametric resonance induced by the onc-photon 
mechanism of EPF creation at the frequency ly ~ 2ra 
(domain R) or as a result of the acceleration of quasi- 
particle e~ pairs in a subcritical external field (do- 
main F). These processes are accompanied by the 
accumulation of EPF density and are limited by the 
condition /(p, t) < 1. For increasing i?o, the domain 
of the momentum space covered by the distribution 
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function is enlarged too ("broadening" of /(p, t), see 
Fig. [7]). This can lead to a rapid growth of the EPF 
density. 

For the strong and moderately strong fields (see 
Fig. [1]) these predictions are based on a realistic KE 
(HI) in the scope of its applicability. For sufficiently 
small fields, Eq <^ Ec, the accumulation R-domain 
comes in the region of the strongly fiuctuating exter- 
nal fields, where the KE (|4]) is not applicable, gen- 
erally speaking, and its solutions have here only ex- 
ploratory character. 

As the EPF density under the accumulation con- 
ditions in the R domain can be very high even for 
relatively weak fields, Eq ^ Ec-, it can be expected 
that the manifestation of the dynamical Schwinger 
effect can be quite possible in this domain due to 
the generation of secondary effects such as, e.g., the 
radiation of annihilation photons. Apparently, these 
remarks can be useful for the discussion of possible 
experiments for observing the dynamical Schwinger 
effect. On the other hand, the accumulation effect 
in the both accumulation domains R and F in the 
short-wavelength domain can be considerable only 
for a sufficiently long duration of the laser field ac- 
tion t^T fFigs. [6l Ffl and fT4|) . This circumstance 
makes difficult an experimental observation of the 
effect in the case of a short single laser pulse. 

In the theory of vacuum EPF creation under the 
action of a periodic laser field of the type ([3]) the 
residual EPF density is the basic final product of the 
theory. However, its definition can be difficult. In 
the present work the residual density is understood, 
as a rule, as the EPF density at the final moment 
of action of some field periods (see, e.g., [13, fisj). 
It corresponds to the density n/p) p3|) . In such a 
case it is implied that the field ^ either prolongs to 
act farther than the periodical field or it switches off 
instantly. Both these cases are rather artificial. It 
would be more satisfactory to consider some pulsed 
field of the type (fT4|) . when a periodical field switches 
on and off gradually or adiabatically. The similar 
situation was discussed in the Sect. IIIII It is very 
important that it does not lead to qualitative change 
in the picture for moderate fields thus demonstrating 
the stability of the obtained results with respect to 
the way of switching the laser field on and off. 

Let us underline also that the approach used in 
the present work is based on a minimal number of 
very general assumptions (spatial homogeneity of 
the external field and its linear polarisation) and is 
in essence an exact consequence of the basic QED 
in the framework of these limitations. This circum- 
stance allows the authors to express the opinion that 
the obtained results in the framework of the given 
approach arc free from some additional approxima- 
tion (e.g., of the WKB type) and therefore a more 
trustworthy in the QC-domain. Moreover, they al- 
low to make first prognoses for the Q-domain. On 
the other hand, the basic kinetic equation (KE) @ 
is valid in the case of a quasiclassical electric field ^ 



only. We extrapolated the domain of its applicabil- 
ity into the Q domain. The foundation for such an 
extrapolation is the assumption that the solution of 
the KE in the boundary domain is valid approx- 
imately because there the quantum fluctuations of 
the electric field are not too large. Thus, we observe 
here the tendency to a non-monotonicity described 
above. In the Q domain the external field can be 
described with the help of the corresponding den- 
sity matrix of the electromagnetic field. Then, the 
quasiclassical field will play the role of a back- 
ground field, i.e. A^^j = -I- A'', where is the 
field operator of the photon component of the total 
field 

For a more adequate investigation of the PV re- 
sponse in the Q domain it is necessary to develop a 
kinetic theory for the generalization of the KE ^ 
taking into account the quantum fluctuations of the 
electric flcld. For understanding the situation, we 
mention the works [45l . l46j in which the vacuum par- 
ticle creation is considered under the influence of a 
stochastic time dependent electric field. We want to 
mention two of the possibilities which emerge here. 
The compressed state of the electromagnetic field 
(e.g., [131 ) is a quantum one but contains a large 
occupation number of photons. As a result, the cri- 
terion (6) of the quasiclassical case can be not ad- 
equate to this situation. For the second example 
we remark that, in the case of a rather strong field 
Eq > 0.1 Ec, the backreaction mechanism becomes 
essential and leads to a stochastic (or close to one) 
internal electric field [1^. Thus, the kinetic descrip- 
tion of vacuum particle creation in the Q domain is 
an actual problem deserving further investigation. 

Considering the muon PV independently of e~e~^ 
PV, it can be expected that the response of the 
PV will repeat on the qualitative level the 
picture described in Sect. IIIII with a shift to the side 
of shorter waves. 

Two regions in the plane of the parameters ([1]) and 
^ remain uninvestigated in the framework of the 
used approach: the top left corner on Fig. [1] where 
short wavelengths and strong fields simultaneously 
persist, and the right boundary of the tunneling do- 
main in the calm valley where A > Aq. It would 
be worthwhile to investigte these areas separately. 
Preliminary results of this work have been reported 
recently [43| and a more elaborate discussion is in 
preparation. 



VI. CONCLUSION 

We have investigated the behavior of the quasipar- 
ticlc EPF generated from the PV under the action 
of a strong laser field. It was shown that particle 
production in the initial stage of the field action is 
characterized by an equilibrium like "thermal" dis- 
tribution. However, at later times the quasiparti- 
clc EPF distribution becomes very complicated and 
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shows a far-from-equilibrium behaviour with distinct 
features depending on the specific domain of the 
landscape. In the subsequent, second part of the 
work we plan to investigate features of the real EPP 
which remains after the laser pulse ceases. More- 
over, the transient phenomenona between quasipar- 
ticle and residual states of the EPP shall be studied. 
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